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ABSTRACT 

The quantitative spectral analysis of medium resolution optical spectra of A 
and B supergiants obtained with DEIMOS and ESI at the Keck Telescopes is used 
to determine a distance modulus of 24.93 ±0.11 mag for the Triangulum Galaxy 
M33. The analysis yields stellar effective temperatures, gravities, interstellar red- 
dening, and extinction, the combination of which provides a distance estimate 
via the Flux- weighted Gravity-Luminosity Relationship (FGLR). This result is 
based on an FGLR calibration that is continually being polished. An average 
reddening of < E(B — V) > ~0.08 mag is found, with a large variation ranging 
from 0.01 to 0.16 mag however, demonstrating the importance of accurate in- 
dividual reddening measurements for stellar distance indicators in galaxies with 
evident signatures of interstellar absorption. The large distance modulus found is 
in good agreement with recent work on eclipsing binaries, planetary nebulae, long 
period variables, RR Lyrae stars, and also with HST observations of Cepheids, 
if reasonable reddening assumptions are made for the Cepheids. Since distances 
based on the tip of the red giant branch (TRGB) method found in the litera- 
ture give conflicting results, we have used HST ACS V- and /-band images of 
outer regions of M 33 to determine a TRGB distance of 24.84 ± 0.10 mag, in basic 
agreement with the FGLR result. We have also determined stellar metallicities 
and discussed the metallicity gradient in the disk of M33. We find metallicity of 
Z & at the center and 0.3 Z & in the outskirts at a distance of one isophotal radius. 
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The average logarithmic metallicity gradient is -0.07 ± 0.01 dex kpc -1 . However, 
there is a large scatter around this average value, very similar to what has been 
found for the H II regions in M33. 

Subject headings: galaxies: distances and redshifts — galaxies: individual (M33) 
— stars: abundances — stars: early type — supergiants 



1. Introduction 



Modern astronomy has advanced to a stage where c qsmological paramete rs are claimed 



to be determined with hitherto unknown precision (e.g. iFrieman et al.ll2008l ). At the same 
time, however, we are confronted with the fact that distances to many nearby galaxies 
are disturbingly uncertain. A n interesting and imp ortant example is the distance to the 
Triangulum Galaxy M33 (see iBonanos et al.l 120061 for a summary). Several independent 
distance determination techniques have been emplo yed for this galaxy since the beginning 



of this decade, includ ing Cepheids and RR Lyrae (IFreedman et al.l 1 199 It iLee et al. 



2007 



Galleti et al 



etary nebulae (ICiarduUo et a 



Saraiedini et al.ll200q). tip of the red giant branch (TRGB) and red c lump stars (jRizzi et al. 



2004 



period variables ( Pierce et al. 



McConnachie et al.ll2004l : iTiede et al.ll2004l: iKim et al 
20041 ). horizontal branch stars (ISaraiedini et al 



200 



2002: 



20021) plan- 
2OO0h . long 



(j), detached eclipsing binary ( Bonanos et al. 2006 ). and 



water masers (IBrunthaler et al.l 120051 ) . The shocking result is that the distance moduli ob- 
tained with these different methods differ by as much as 0.6 magnitudes, which is more than 
30 percent in linear distance. Even when the water maser distance were discounted in this 
comparison as justified by its large intrinsic uncertainty of 0.45 mag, the spread in distance 
moduli is still about 0.4 magnitudes, significantly higher than the uncertainties quoted in 
the individual works. 

A significant part of the uncertainty comes from interstellar extinction. M 33 is a mild- 
inclination star-forming galaxy with a patchy dusty disk that exhibits inhomogeneous intrin- 
sic extinction. This has a considerable effect for many of the aforementioned methods. In 
particular, distance determination using Cepheids at V and / bands is severely affected, as 
they a re part of a young stell ar population in a dust y star- forming environment. As pointed 
out by lKim et al.l (120021 ) and ICiarduUo et al.l (120041 ). the difference between the shorter dis- 
tances obtained with Cepheids and the longer distances found with other methods would 
disappear given s maller reddening than assumed in the Cepheid work cited above. Interest- 
ingly, the work by lBonanos et al.l (120061 ). which uses quantitative spectroscopy of a detached 
eclipsing binary and, thus, comes up with an accurate determination of reddening and ex- 
tinction, yields a large distance modulus of 24.92 mag, 0.3 magnitudes larger than the one 
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obtained with V- and /-band photometry of Cepheids. 



The solution to resolving this discrepancy is to use a method that is either unaf fected by 
redde ning and extinction uncertainties such as K-band photometry of Cepheids (e.g. iGieren et al. 



20051 ). or one that includes the accurate determination of reddening for each individual star 
gauged in the determination of the distance modulus. An example for the latter involves the 
quantitative spectroscopy of massive blue supergiant stars and the use of the Flux- weighted 
Gravity-Luminosity Relationship (FGLR; lKudritzki et al.ll2003l ) . In this paper we will apply 
this method to determine the reddening, extinction, and distance to M33. 

Blue supergiants of B and A spectral types are the optically brightest stars in the uni- 
verse with visual magnitudes up to My = —9.5 mag. Their quantitative spectral analysis 
based on state-of-the-art non-LTE radiative transfer models allows for precise determina- 
tion of stellar parameters, metallicity, re ddening, and extinction (jUrbaneja et al.l l2005al ; 
Przybilla et al.ll2006l ; iKudritzki et al.ll2008l ). Then, with stellar temperatures T e g and gravi- 
ties log g, one can use the de-reddened apparent stellar magnitudes and the tight relationship 
between flux-weighted gravity log gp = log g — 4 log (T cS x 10~ 4 ) and absolute bolometric 
magnitude M^ Q \, to determine a distance via the FGLR. 



The concept of FGLR has been introduced by IKudritzki et al.l (120031 ). It is based on 
the assumption that massive stars evolve through the Hertzsprung-Russell (H-R) diagram 
temperature domain of B and A spectral types at constant luminosity and mass. In a 
detailed spectroscopic stud y of blue supergiants in the nearby spiral galaxy NGC 300, which 
is at a distance of 1.9 Mpc, IKudritzki et al.l (120081 . hereafter K08) v erified this concept and 
contributed towards the calibration of the FGLR. lUrbaneja et al.l (120081 ) then applied the 
FGLR method and determined the distance to the metal-poor, irregular Local Group galaxy 
WLM, the result from which agreed well with most recent distance determinations based on 
K-band photometry of Cepheids and I-band photometry of the TRGB. 

Motivated by the success of previous application of the FGLR as well as the prospects 
of this work, we undertake an independent investigation of the distance to M33 by means of 
a quantitative analysis of medium resolution optical spectra of blue supergiant stars. 

This paper is organized as follows: In §2] we describe the data sample used to obtain 
the FGLR of M33 as well as the observations undertaken at Keck Observatory. In £J3] we 
detail the process of the quantitative spectroscopic analysis. We present the results of the 
spectral analysis in §H and discuss stellar parameters and stellar evolution, metallicity and 
metallity gradient, and reddening and extinction. This provides the basis for a new distance 
determination using the FGLR in £j5j A discussion of the new result compared with distances 
determined from other methods is given in §5.31 Finally, we present our conclusions and 
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future outlook in £|6j 



Observations 



Distance determi nation using the FGL R method requires a sufficiently large sample (10 
to 20 objects, see K08. IUrbaneja et al. boost ) of supergiants with good multi-band photometry 



and spectra of reasonable quality that are suitable for quantitative analysis using model 
atmosphere and radiative transfer techniques. We pre-selected a sam ple of such potentia l 



targets that fulfill the initial color-magnitude criteria described in K08. iMassey et al.l (120061 ) 
provides an excellent source for photometry of all these objects; spectra were collected with 
the Keck Telescopes on Mauna Kea. In our first observing run for this project on November 
1, 2003 with good seeing (0'/6), we used the DEIMOS spectrograph with a 1200 1/mm grating 
yielding a dispersion of 0.33 A/px and a spectral resolution (FWHM) of 1.6 A. The total 
integration time was 135 and 93 minutes for the blue and red spectra, respectively. A total 
number of 49 targets were observed with one DEIMOS mask. Among these, we selected 10 
objects with spectral type ranging from B3 to A3 based on the requirement of a minimum 
signal-to-noise ratio (SNR) needed to carry out the ensuing spectral analysis. 

Additional observations of 7 more A-type supergiants were carried out under good 
weather conditions (clear sky, ty pical seeing around l'.'O) on October 8, 2005 and October 4, 



2007 with the ESI spectrograph (ISheinis et al.ll2002l ). The instrument was used in echellette 



mode, with a slit width of 0"75, providing R~9000 (equivalent to 0.6 A FWHM at 5500 
A). Depending on the visual magnitude of the object, one to three exposures were acquired, 
with the exposure time adjusted to provide a final SNR ~ 100 in the continuum. Unlike in 
the case of DEIMOS, the ESI stars were selected to cover a wide range of galactocentric 
distances. 

To further increase the sample size for FG LR analysis and calibr ation, we included 



additional stars from the spectroscopic study by lUrbaneja et al.l (l2005bl ). In this work, 11 
B supergiants of early spectral type (B0 to B2) were analyzed quantitatively using state- 
of-the-art non-LTE model atmospheres to investigate the M33 metallicity gradient. While 
these stars were well suited for the abundance study, however, care must be taken while 
considering their selection as targets for the FGLR distance determination method. Two 
of them (B133 and 1054) were found to be multiple systems with companions that would 
affect magnitudes and colors. Two others (110A and OB10-10) turned out to be extreme 
blue hypergiants, with very high luminosities, rendering the determination of their surface 
gravities less reliable. For object 1137, only an upper limit for the mass loss rate could be 
found, which means that, because of wind contamination in the Balmer lines, only an upper 



- 5 - 



limit for the gravity could be determined. The spectra for object 0900 suffered from strong 
H II region emission and consequently-distorted profiles of the two observed Balmer lines H7 
and B.5, again leading to a rather uncertain stellar gravity. The remaining five objects have 
very well determined stellar parameters and have been incorporated as targets in the FGLR 
distance determination technique. This brings the total number of M33 supergiants in our 
sample up to 22. The list of all the objects is given in Tbl. [TJ 



The new ESI data were reduced in the same way as described in lUrbaneja et al.l (l2005bl ). 
On the other hand, the DEIMOS spectra were initially reduced using the DEEP2 Reduc 

tr 



tion Software (e.g. iMarinoni et al.ll200ll ). but the resulting SNR in the short wavelength 
range were insufficiently low because the subroutines in the pipeline optimizes in the longer 
wavelength regime. Therefore, the data reduction of the DEIMOS data was carefully redone 
step-by-step using the standard Image Reduction and Analysis Facility (IRAF) packages. As 
demonstrated in Fig. HJ the final extracted spectra have a clearly improved SNR in the blue 
region, giving better results than the DEEP2 pipeline did. The noticeably-improved SNR 
allows us to identify and measure the relevant spectral features more precisely. However, 
even after this improvement in SNR, the poor blue sensitivity of DEIMOS (as can be seen 
in Fig. [T]) turned out to be an extra challenge during the spectral analysis. 



3. Spectral Analysis 



The analysis method to be applied to the individual targets depends on the targets' spec- 
tral types. For early B supergiants, we use the Si il/lll/iv equilibrium based on metal-line- 
blanketed non-LTE model atmospheres and line formation calculations that i nclude effects of 



stellar winds and spherical extension to constrain effective temperatures (see lUrbaneja et al. 
2005aB ). For late B and early A spectral types, the determination of effective temperature 



utilizes the Balmer jump or ionization equilibrium such as Mg 1/11 and O 1/11. These mod- 
els involve LTE line-blanketed model atmospheres in hydrostatic equilibri um with plane- 



paral l el geometry and very detailed non-LTE line formation calculations (IPrzybilla et al. 



20061 ; iKudritzki et al.l 120081 ). In consequence, as a first step towards quantitative spectral 
analysis, we determined spectral types for all 14 newly-observed objects through a qualitative 
comparison with a complete set o f Milky Way and SMC supergiant template spectra (see 
Bresolin et al.ll200ll ; iLennonl 119971 ). Only three objects were found to exhibit early spectral 
type of B3 while that of the rest ranges between B8 to A2 (see Tbl. |2J). 



Spectra analyzed in the FGLR studies by K08 and lUrbaneja et al.l (120081 ) had a rather 
low resolution of 5 A. While this resolution is sufficiently good for determining stellar grav- 
ities from the Balmer lines, it precludes the use of ionization equilibria for temperature 



- 6- 



determination in the case of late B and early A supergiants because the weak spectral lines 
of the neutral species disappear in the noise of the spectra. As an alternative, the Balmer 
jump or the strengths of the neutral helium lines were utilized instead. However, with the 
spectroscopic resolution in our Keck DEIMOS and ESI spectra (1.6 a nd 0.6 A respectively) , 
we can employ the technique of ionization equilibria in this work (see iPrzybilla et al.ll2006l ). 



The basis for the spectral analysis of late B and A type supergiant stars is the grid of 
model atmospheres and line formation calculations in K08 that spans an extensive area in 
(T eff , log g) parameter space (see Fig. 1 in K08) and covers a large range in metallicity from 
0.05 Zq to 2 Z & . We refer to K08 for a description of the computed grid (abundances, choice 
of microturbulence velocities, etc.). T he input physic s of th e model atmospheres and line 
formation calculations are detailed in IPrzybilla et al.l (120061 ) . Regarding the mi d B-types, 



we u se the unified non-LTE model atmosphere/line formation code fastwind (IPuls et al. 



20051 ) to create a small grid of models, with the parameters properly selected to explore the 
expected small parameter space covered by these three objects. 

The analysis proceeds in three steps. First, projected rotational velocity is determined; 
second, stellar effective temperature and gravity; third, stellar metallicity. If needed, steps 
two and three are iterated. In the following, we describe the individual steps. 



3.1. Rotational Velocities 



Projected rotatio nal velocities v sin i of late B and early A supergiants are usually be- 
tween 30 to 70 kms -1 (jHowarthl|2004j ) comparable to or somewhat larger than the resolution 
of the DEIMOS and ESI spectrographs, respectively. In other words, we have to determine 
the amount of rotational line broadening in addition to instrumental broadening for each 
individual star before we start the spectral analysis via comparing synthetic with observed 
line profiles. 

With the large number of objects to be analyzed, it is important to have a simple 
method to accurately estimate v sin % that is independent of stellar parameters as well as the 
details of the final atmospheric model synthesizing the observed spectrum. A good approach 
is to look at lines that are weak enough to not suffer from pressure broadening to a large 
degree but s till sufficient ly strong such that the rotation-broadened profile shape can be well 
determined JGravi Il992h . The Mg II A4481 line, the Si II A6347, 6371 lines and the Fe II 
A4515, 4535 lines are ideal candidates for late B through early A supergiants while the Si III 
lines at A4552, 4567, 4574 are useful for early B supergiants. For those lines we construct 
a Gaussian profile with very narrow FWHM (0.1 A) that has the same equivalent width as 
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the observed line. We then convolve this initial profile with the instrumental profile and a 
rotational profile, where we steadily increase v sin i. An example in the case of the DEIMOS 
spectra is given in Fig. [2J A x 2 minimization of the difference between the observed and 
calculated profiles in turn indicates projected rotational velocities. We note a caveat that 
Mg II A4481 is a blend of two strong components separated by 0.2 A; determination of wsin % 
with a single Gaussian may tend to overestimate v sin i in these cases. 

The lower spectral resolution and SNR achieved for the DEIMOS data do not allow 
us to disentangle rotation and macroturbulence. Therefore, the value of v sin i thus ob- 
tained incorporates macroturbulence and might overestimate the true rotational velocity. 
This, however, will not affect the determination of stellar parameters and subsequent FGLR 
distance analysis. The situation is more optimistic for the late B and early A supergiants 
observed with ESI that have significantly better spectral resolution and SNR; we thus de- 
cided to apply a two-step procedure. First, v sini was determined by a Fourie r trans form of 
individual lines (such as the ones mentioned previously) as described in I Gray! (119921 ). With 
this value fixed, the macroturbulence velocity was th e n esta blished by a by-eye fit of the ob- 
serve d profile similar to that done in iPrzybilla et al.l (120061 ) . The early B supergiant sample 
from lUrbaneja et al.l (j2005bl ) were in a situation similar to that of our DEIMOS spectra: in 
order to obtain a reasonable SNR, the ESI spectra had to be binned to the detriment of the 
resolution, which did not allow then to discern between rotation and macroturbulence. 



A list of rotational and macroturbulence velocities for our sample of supergiants is given 
in Tbl. [2J The uncertainties are on the order of 5 kms -1 . 



3.2. Gravity and Effective Temperature 

The standard technique for determining gr avities and effective te mperatures of late B 



and early A supergiants is described in detail by IPrzybilla et al.l (120061 ) and K08. We apply 



a very similar but somewhat more efficient technique in this work. 

The first step involves using information about gravity and effective temperature that 
are provided by a fit of the Balmer line profiles. Balmer lines are an excellent indicator of 
stellar gravity because of their strong density dependence on pressure broadening through the 
Stark effect. An example is given in Fig. 0for object No. 11 of our sample. At a fixed effective 
temperature, the step size of each gravity fit on the grid is 0.05 dex. Since the excitation 
of the second level of hydrogen and the ionization of hydrogen is temperature-dependent, 
the strength of the Balmer lines also depends on effective temperature. This means that the 
Balmer lines may be fitted just as well as higher (lower) effective temperature and higher 
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(lower) gravity than shown for the example in Fig. [3j The corresponding fit curve on the 
(log g, T c g)-plane is demonstrated in Fig. HI Since the major motivation for this work is 
a determination of the distance to M33 using the stellar flux-weighted gravity loggp, we 
also plot the same fit curve for the Balmer lines on the (log gp, T eff )-plane. We note that 
for stars sufficiently hot, log gF can be solely determined by fitting the Balmer lines. The 
physics behind this phenomenon is discussed in detail in K08. 

While the analysis of the Balmer lines provides a strong constraint for gravity and 
temperature, we need a second, more temperature-dependent spectral feature to further 
constrain the effective temperature along the Balmer line fit curve in Fig. |H The most 
accurate s olution is to use ioniza tion equilibria such as Mg i/ll for late B an d early A super- 



giants (see iPrzybilla et al.ll2006l ) and Si ii/m/iv for early B supergiants (see lUrbaneja et al. 



2005bl ). Fig. demonstrates how this ionization equilibrium changes with (T e g, log g) pairs 



and sets the best-fit effective temperature of object No. 5 to be 8750 ± 250 K. 

Because of the relatively high resolution and high SNR (~100), the ionization equi- 
librium method works extremely well for all our targets with ESI spectra. The DEIMOS 
spectra exhibit lower resolution and lower SNR (~ 30-70), but the method is still successful 
for many of the remaining targets, as shown in Fig. [5j However, there were some DEIMOS 
targets for which no Mg I lines could be observed because either the stars were too hot 
or the spectra were too noisy, or both. For such cases we applied an alternative technique 
introduced by K08 that made use of the very strong temperature dependence of the neu- 
tral helium lines in the temperature range from 9000K to 13000K. An example for such a 
DEIMOS case is given in Fig. [6j where a fit of the He I lines of object No. 1 yields T c g = 
10000 ± 500K. The fit of the Balmer lines for the same object is also presented in the figure. 
Naturally, using the He I lines requires making an assumption about the helium abundance, 
which is slightly enhanced in the atmospheres of late B and early A supergiants but varies 
only between y = N(He)/(N(He) + N(H) ) = 0.11 to 0.1 3 acco rdin g to the high resolu 



tion, high SNR quantitative spectroscopy by IPrzybilla et al.l (120061 ) and ISchiller &: Przybilla 



( 120081 ). We have adopted y = 0.12 for the fit in Fig. [61 We note, though, that the influence 
of temperature variation on the Hel lines in this temperature range is much stronger than, 
for instance, changing y from the solar value of 0.09 to 0.13. We also note that for all cases in 
this temperature range, where we were able to constrain the effective temperature with ion- 
ization equilibria (e.g. all our objects observed with ESI), the agreement of the observed He I 
lines and computed models is excellent, lending support to the reliability of He I method. 

Tbl. [2] gives the effective temperatures and gravities thus obtained along with their 
associated uncertainties. Note that in most cases the errors for log g are larger than that 
for log gp-The reason is that the gravity fit curve has a stronger dependence on temperature 
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than the flux-weighted gravity fit curve, as illustrated in the upper and lower panels of Fig.HJ 
respectively (see K08 for a detailed discussion of this effect). 



3.3. Metallicity 

The use of ionization equilibrium such as Mg i/ii as illustrated in Fig. presumes 
a stellar magnesium abundance or, at least, a stellar metallicity (the dependence of the 
Balmer or Helium lines on metallicity is very weak and can be neglected). As a starting 
point for T e g and log g determination, we have a dopted metallici t ies in a greement with the 



metallicities and the metallicity gradient found by lUrbaneja et al.l (l2005bl ). After constraints 
have been placed on the values of T e g and log g, we then check whether the initial metallicity 
assumption was correct. This process is exemplified in Fig. [7] and detailed as follows. We 
define spectral windows, which are dominated by metal lines and for which a good definition 
of the continuum is possible. We then select all available models with different metallicities 
(typically ranging from \Z\ = -1.3 to 0.3 dex, with [Z] = log Z/Z Q ) at the corresponding (T eff , 
log g) and, for each of these models, carry out a pixel-by-pixel comparison of the calculated 
and normalized fluxes within each spectral window. The metallicity over each spectral 
range is thus the metallicity associated with the minimum % 2 -value; the average [Z] from all 
windows then gives an estimate of the stellar metallicity. If the value of [Z] thus obtained 
differs from that used for the ionization equilibrium by more than 0.1 dex, we reiterate the 
T e g and log g determination process using this new [Z\ until the metallicity measurements 
converge. This procedure is repeated for models at the extremes of the error box for T e g and 
log g, which, together with the dispersion in [Z\ obtained in all spectral windows, defines 
the uncertainty in [Z] (see K08 for details). The stellar metallicities with uncertainties for 
our sample of supergiants are given in Tbl. [21 Note that the metallicities obtained via this 
technique for the late B and early A supergiants reflect mostly the abundances of heavy 
elements su ch as iron and t i tanium . For the early B supergiants, the same method as was 
described in lUrbaneja et al. (h()05bl ) was applied: the metallicities are an average of oxygen, 



magnesium and silicon abundances. 



3.4. Interstellar Reddening, De-reddened Magnitudes, Stellar Luminosities, 

Radii, and Masses 



With effective temperature, gravity, and metallicity constrained from our spectral anal- 
ysis, we know the intrinsic energy distribution of our target stars from the flux distribution 
of the model atmosphere computed with the final stellar parameters (see iPrzybilla et al. 
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2006 ; iKudritzki et al.ll2008l ; lUrbaneja et al.ll2008l ). Comparing with the obser ved fluxes from 



Massev et all (120061 ) and applying the interstellar extinction law found by ICardelli et al. 



( 119891 ). we can obtain an accurate estimate of interstellar reddening E(B — V) for each 
individual supergiant (see K08 for details). Assuming Ay = 3.1 E(B — V) for the relation- 
ship between reddening and visual extinction yields the de-reddened apparent bolometric 
magnitude m^ \ given by 

iTiboi = m v — A v + BC , (1) 

where my is the observed apparent visual magnitude; BC is the bolometric correction, 
which is also given by the flux distribution of the model atmosphere calculated with the final 
stellar parameters. Note that K08 give an accurate analytical formula to compute BC as a 
function of T e g, log g, and [Z]. Stellar photometry including intrinsic colors, reddening, and 
bolometric correction is summarized in Tbl. [3j 

Assuming an appropriate distance modulus, for instance as determined via the FGLR 
method (see below), the de- reddened apparent bolometric magnitude mt, i can subsequently 
be used to evaluate the absolute bolometric magnitude Mboi or, equivalently, the stellar 
luminosity log (L/L & ). Stellar luminosities coupled with effective temperatures then yield 
stellar radii. Tbl. H] lists the absolute bolometric magnitudes and radii obtained in this way. 

As explained by K08, there are two ways to determine stellar masses. We can use 
stellar gravities together with the radii to directly calculate masses; masses subsequently 
determined are referred to as spectroscopic masses. Alternatively, masses can be esti- 
mated by comparing the locatio n of our target stars on the H-R diagram with evolutionary 
tracks (IMaeder &: Meynetll2005l ); these masses are called evolutionary masses. Both masses 
along with their uncertainties are given in Tbl. HI 



4. Results 

In this section we discuss the main results of our quantitative spectral analysis that are 
compiled in Tbl. Tbl. El and Tbl. ffl 



4.1. Interstellar Reddening and Extinction 

Fig. [8] shows the distribution of interstellar reddening E{B — V) among the stars in our 
sample. We find a wide range from 0.01 to 0.16 mag with an average value of < E(B — 
V) > = 0.083 mag. The individual values are significantly lower than E(B — V) = 0.20 mag, 
the reddening value as adopted in the HST distance scale Key Project study of Cepheids 
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by iFreedman et al.l (ll99ll ) (see also lLee et al.ll2002l ). This will have important repercussions 
for the discussion of the distance to M 33 in §5.31 We note that our average value is higher 
than the foreground value of E (B — V) = 0.04 mag, which has been derived from the reddening 
maps by ISchlegel et al.l (119981 ) . We also note that the wide range in reddening found in this 
study is similar to the range in NGC300 derived by K08 and demonstrates the need for 
reliable individual reddening determinations. 

An independent method to acquire information about interstellar reddening in M33 is 
the study of the Balmer decrement of H II regions. While H II regions and the associated stars 
within are generally younger than A and B supergiants as well as Cepheids, which in many 
cases have already migrated into the field, and while the average reddening of H II regions is 
hence very likely somewhat hig her, it is still useful to d i scuss E(B — V) values obtained with 
this technique. Most recently, iRosolowsky fc Simonl (120081 ) in their comprehensive study 
of M33 H II regions have published reddening values c^Hp) for a large sample. Adopting 



E(B — V) = 0.676 c(Hfl) (iMagrini et al.l l2007bl ). we produce Fig. [9] with reddening values 
from IRosolowsky fc Simonl (I2008T). where we plot E (B — V) as a function of angular galactic 



distance, in units of R25 ( Ide Vaucouleurs et al.lll995l ). The result is very informative: there is 
a large scatter in reddening at all distances; however, the average E(B — V) value (excluding 
the three most extreme cases with E(B — V) >0.6 mag) is < E(B — V) > = 0.11 mag, 
which is in close agree ment to what has been found for our blue supergiants. We note that 
Magrini et al.l (j2007bl ) quote an average reddening of < E{B — V) > = 0.22 mag from their 
study of H II regions; however, their work systematically overestimates the reddening value, 
since, when measuring Balmer emission line fluxes and the Balmer decrement, they have not 
corrected for the effect of the underlying stellar Balmer absorption. 



4.2. Stellar Properties 
Fig. [10] (upper panel) sho ws the location of all t argets on the (log g, log T e ff)-plane 



along with evolutionary tracks (IMeynet fc Maederll2003l ) calculated for solar metallicity that 



include effects of rotation and anisotropic mass loss. In this diagram, the data points of 
the stars stand independent of any assumption about the distance to the galaxy and rely 
completely on quantitative stellar spectroscopy. There is an indication of two evolutionary 
sequences. The six early B supergiants (No. 17 to 22 on our target list) and one late 
B supergiant (No. 16) align with the 40 M & evolutionary track, while the others appear 
to be in the 20 to 25 M & range. One object (No. 4) that coincides with the 15 M Q 
evolutionary track seems to be of lower mass. The fact that all the early B supergiants are 
more massive is a selection effect: the objects have been chosen according to their apparent 
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visual magnitudes. Since the early B supergiants are hotter than the rest of the sample, 
their bolometric corrections and hence their luminosities and, consequently, their masses are 
larger (see discussion in K08). 

Complementary information about the evolutionary status of the targets can be derived 
from their location on the H-R diagram. With a distance modulus of \x = 24.93 ± 0.11 mag 
as determined in $5j we have calculated stellar radii and luminosities as described in the 
previous section. This allows us to position all targets on the H-R diagram in another 
comparison with the same evolutionary tracks (see Fig. [TUJ bottom panel). The advantage 
of using this diagram is that the theoretical dependence of stellar luminosity on mass is very 
strong and, given the size of the error bars, allows for a more accurate assessment of the 
original stellar masses. Most obviously, the same objects that appear to be more massive in 
the previous plot are also the more luminous ones associated with higher-mass evolutionary 
tracks. This nicely confirms the purely spectroscopic analysis, at least qualitatively. 

However, we encounter a problem at lower luminosities. Now there are three objects 
(instead of just one) located on the 15 M Q track at lower effective temperature. They are 
No. 2, 4, and 5 in order from lower to higher luminosity. Target No. 2 is particularly 
suspicious: it has the lowest luminosity of all targets indicating an original mass of less than 
15 M , but on the (log g, log T e ff)-plane it has the lowest gravity and effective temperature, 
resembling an object of significantly higher mass. This indicates a discrepancy between the 
spectroscopic result and stellar evolutionary path for this specific object. 

The discrepancy becomes more apparent when we plot spectroscopic stellar masses 
against stellar luminosities and compare the results with the corresponding relationship 
derived from currently available stellar evolutionary tracks at the effective temperatures of 
early B supergiants and A supergiants (see K08 for a simple fit formula). This is done in 
Fig. [HJ which reveals that the spectroscopic mass of target No. 2 (the object with the lowest 
mass and luminosity) does not fit on the stellar evolutionary mass-luminosity relationship. 
Its mass is clearly too low for its present luminosity. We note that K08 in their study of blue 
supergiants in NGC300 have found a very similar outlier object (their target No. 17). They 
speculated that these objects could be stars of lower initial mass now evolving back from 
the red giant branch to become blue sup ergiants again at s i gnifica ntly higher luminosity. 



However, an inspection of the tracks by iMeynet fc Maederl (120031 ) does not support this 



speculation. None of the tracks with 12, 15, and 20 M & loops back from the red giant 
branch, whereas the tracks with significantly lower masses do not reach the luminosity of 
object No. 2 within their paths. A more attractive explanation appears to be close binary 
evolution during which No. 2 (and No. 17 of K08) might have lost a significant fraction 
of their original mass while roughly keeping their original luminosity. Expectedly, such a 
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scenario will need to be confirmed by detailed evolutionary calculations. Regardless of the 
details, it seems that this object has not followed the same evolutionary path as the rest 
of our sample. It is obvious that the star will not fit into the FGLR defined by regular 
supergiant stars. Therefore, for distance determination using the FGLR-method we will 
exclude object No. 2, in the same way K08 excluded their object No. 17. 

Fig. [11] seems to indicate that the observed relationship between spectroscopic masses 
and luminosities for the remaining objects is in slight disagreement with stellar evolution 
models since most of the objects are somewhat overluminous for their mass. In order to 
discuss this more quantitatively we use the stellar evolutionary mass-luminosity relationships 
for A and B supergiants in Fig. [11] (see also the fit formulae in K08) to convert observed stellar 
luminosities into evolutionary masses (as described in §3]). In Fig. [T2]we plot the logarithmic 
ratio of spectroscopic to evolutionary masses as a function of luminosity. As we can see 
from these plots and Tbl. H] spectroscopic masses are on average 0.06 dex lower than their 
evolutionary counterparts (not including No. 2 in calculating the average), indicating indeed 
a small systematic discrepancy between spectroscopic results and stellar evolution theory. We 
note that while evolutionary masses have much smaller error bars than spectroscopic masses 
do due to the steep dependence of luminosity on mass, they might be subjected to systematic 
uncertainties such as the influence of mass loss over the stellar lifetime or rot at ionally- induced 
mixing effects on theoretical luminosities. Of course, spectroscopic masses might be subjected 
to systematic uncertainties as well, in particular through the way of pressure broadening of 
the Balmer lines is being accounted for in the process of spectroscopically determining stellar 
gravities as demonstrated in the previous sections. 

In order to assess whether or not the discrepancy encountered is specific to our M33 
sample, we repeat the logarithmic ratio of the spectroscopic to evolutionary masses analysis 
for stars in other galaxies in the bottom panel of Fig. O This time we examine targets from 
the eight galaxies used fo r the F GLR calibration by K08 and from the most recent study of 



WLM by lUrbaneja et al.l (120081 ). We do not see any indication of a special systematic effect 
that is only related to M33. The average logarithmic ratio for the 66 targets overplotted is 
-0.04 dex. 

Due to the fact that both evolutionary and spectroscopic masses depend on a presumed 
luminosity, it is worthwhile to discuss whether or not an inappropriate choice of the distance 
modulus would cause a systematic effect. Since M spcc oc gR 2 oc d 2 , L oc d 2 , and a mass- 
luminosity relationship of the form L cx M® vol , we derive a logarithmic ratio of spectroscopic 
to evolutionary masses with dependence on a and d as the following: 

log M spcc /M cvol = (2 - 2/ a) log d . (2) 
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With a ~ 3.5 (see Fig. [TTJor K08), we obtain logM spec /M evo i = 1.4 log d, which means 
that the distance modulus for M33 would have to be 0.2 mag larger than what we have 
adopted. As the discussion of the distance to M33 in §5] will show, this is unlikely. 



4.3. Metallicity and Metallicity Gradient 



Findings from our quantitative spectroscopic study can be used to assess the metallicity 
of the young stellar populations in M33 as well as the results of emission line studies of 
H II regions. Fig. [13] displays the stellar metallicities obtained in consequence as a function 
of dimensionless angular galactocentric distance. In this, and the corresponding following 
plots, we use i? 2 5 = 35.40' (jde Vaucouleurs et al.lll995l ). Our basic outcome is that the stellar 
metallicity is close to the solar value at the center and that it decreases by a factor of three 
at the galactocentric distance of an isophotal radius. A simple linear regression of the form 



[Z] = [Z] + [Z^iR/R 



25 j 



(3) 



0.73 ±0.09 dex i?^ 1 for 



yields [Z] = 0.09 ± 0.04 dex for the central metallicity and \Z\\ 
the angular gradient (equivalent to -0.07 ± 0.01 dexkpc -1 for the distance used in this paper). 
We note, however, that the scatter is significant around 0.5i?2s; it is very likely a real 
phenomenon rather than an artifact of the uncertainties of our analysis. 

Fig. dH is similar to F i g. [TBI but it also includes the H II region Oxygen abundances 
compiled by IMagrini et al.l (|2007al ) , H II Neon abundances as determined mos t recen tly by 
Rubin et al.l (120081 ) . and Cepheid metallicities as obtained by iBeaulieu et al.l (120061) . (We 
found that the de-projected galactocentric distances given by IBeaulieu et al.l (120061 ) contained 
errors that we have corrected for Fig. [T4"l) . The value for the solar Oxygen abundance adopted 
to normalize the H II dat a bv IMagrini et all J2007ah is log [N(0)/N(H)) +12 = 8.69 dex 
( lAUende Prieto et al.ll200ll ). From this sample, we considered only those H II regions for 
which a determination of the electron temperature was possible. Since the value of the 



solar Neon abundance seems to be rather uncertain (see discussion in iRubin et al.ll2008l ). 
we normalize Neon to the average yalue for B stars in the solar neighborhood and adopt 
log [N(Ne)/N(H)} +12 = 8.08 dex JPrzvbilla et alihoosh . 



The different data sets seem to converge: they all indicate a significant abundance 
gradient. Combining both sets of H II region abundances, we derive a slope of -0.55 ±0.09 
dexi?^ 1 (-0.06 ± 0.01 dexkpc^ 1 ), statistically consistent with the stellar results. With respect 
to the intercept, a comparison is less robust since several issues could be affecting the zero 
point of the different abundance scales. The nebular data provide -0.11 ±0.04 dex, slightly 
below the stellar result, but nevertheless consistent and not unexpected given the systematics 
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that could be affecting both samples (stars and ionized gas). We regard these as consistent 
results. A least-square fit to all objects in this diagram yields [Z] Q =-0.01 ±0.03 dex for 
the central value and [Z]i =-0.71 ± 0.07 dexi?^ 1 for the gradient. Within the uncertainties, 
this agrees with the results obtained from the supergiant stars alone. We note again that 
the results for the H II regions displayed in this figure also indicate a large scatter around 
the regression curve. 

This significant metallicity gradient is in obvious disagreement with the result found 
by iRosolowsky fc Simon! (120081 ). who carried out a comprehensive study of H II regions in 
their M33 metallicity project. Fig. [151 compares the Oxygen abundances of our supergiants 
with that of their H II r egions. While our data poin ts largely coincide with theirs at large 
galactocentric distances, IRosolowsky fc Simon! (120081 ) find many H II regions with very low 
oxygen abundance in the central region of M33 for which we have no stellar counterparts in 
our sample. At this point, we have no explanation for this discrepancy. 



5. Distance 

In this section we employ two different methods to determine the distance to M33. First, 
we use the stellar parameters from this spectroscopic study and apply the FGLR method; 
second, we analyse Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) 
photometry of three outer fields of M 33 and obtain a distance from the J-band magnitude 
of the tip of the red giant branch (TRGB method). The results from each technique will be 
discussed accordingly. 



5.1. FGLR and Distance to M33 



The FGLR is a tight correlation between the flux-weighted gravity (gp = g/T 4 e g) 
and the absolute magnitude of BA supergiants. The ph ysical background, det ection and 
the calibration in nea rby galaxies has been described in iKudritzki et all (120031 ) and K08. 
Urbaneja et all (120081 ) were the first to use the FGLR for distance determination and found 
a distance modulus of 24.99 mag for the metal-poor dwarf galaxy WLM, in good agreement 
with most rece nt TRGB distance de terminations. Here, we follow the same procedure as 
was detailed in lUrbaneja et all (120081 ). 



The FGLR has the form 



M bol = a(log g F - 1.5) + b 



(4) 
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with the most recent calibration provided by K08, a = 3.41 and b = -8.02. Our spectroscopic 
analysis provides de-reddened apparent bolometric magnitude m^ \ (see calculation details in 
§3.4) and flux-weighted gravity for each of our targets (except for object No. 2, as explained 
in §4.2) that we can fit with a regression of the form 

^bol = flM33(log QF - 1-5) + &M33 • (5) 

The fit result is shown in Fig. [16] (top). Since our M33 targets span only a limited range 
in g F compared to the K08 sample, we adopted the slope value provided by K08 (letting 
«M33 = a ) and fit only &M33- The difference between b and &M33 yields the distance modulus, 
which we determine to be /1 = 24.93 ±0.11 mag (the error is calculated similarly as in 



Urbaneja et al.ll2008l ). Fig. [16] (bottom) then compiles the FGLR results for our M33 targets 



using this distance modulus as well as that for 9 other galaxies from K08 and lUrbaneja et al. 



(120081 ). The plot does not indicate any systematic differences between the M33 sample and 



the remaining sample. 

We note that the present calibration of the FGLR by K08 is based on a sample of 
supergiants selected from NGC 300 and seven Local Group galaxies for which distances 
were adopted. A cleaner, more accurate and systematic method to calibrate the relationship 
would be to use a large sample of blue supergiants in the Large Magellanic Clouds (LMC). 
This work is presently under way. 



5.2. TRGB and Distance to M33 



We use observations of M 33 available through the HST archive to measure a distance 
based on photometry of stars at the TRGB. The footprints of the observations taken with 
ACS are shown in Fig. [17] (top) as yellow outlines. The two fields towards the east (left) side 
of the image are from HST Program 9479's study of M33's halo; these fields are combined 
in the ensuing analysis. The remaining field to the south of the image comes from HST 
Program 10190's observations of the outer disk. We choose these sets of images from among 
the many available through the HST archive in order to balance the need of minimizing issues 
of stellar crowding and reddening within M33 with the need for a well-populated RGB. 

Both programs use the F814W and F606W filters, which can be described as 'wide 
V and 'wide V filters, respectively. We perform photometry and conduct artificial star 
tests on t he images usin g the DOLPHOT software package, which is a modified version of 
HSTphot (lDolphinll2000l ). The resulting color- magnitude diagrams are presented in Fig. [T7J 
(bottom). Using the real and artificial photometry we measure the - F814H^ magnitude an d 
color of the TRGB using a maximum-likelihood method described by lMakarov et al.l (120061 ). 
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The F814W magnitude of the TRGB measured in the halo fields from Program 9479 
is 20.92+°;°? with FQOQW - F8UW = 1.35 ± 0.05, while the observations of the disk from 
Program 10190 yields 20.89±g;g§ with F606W - F8\m = 1.47±g;g|; the TRGB magnitudes 
are marked with broken horizontal lines in Fig. [TT1 The flux from stars at the TRGB 
is least sensiti ve to age and m et allicity in t h e /-b and, and has an absolute magnitude of 
Mj ~ —4.05 (IRizzi et all 120071 ). iRizzi et al.l (120071 ) also describe a method for calculating 
a distance modulus that includes a zero-point calibration and metallicity correction; the 
formula for which we reproduce below: 



DM = m TRGB - A F8UW + 4.06 - 0.20[(F606W - F8UW) - (A F606W - A F8UW ) - 1.23] , (6) 



where w-trgb is the apparent magnitude in F814W, and A F81iW and A F q 06W are foreground 
extinction fo r the two filters, resp ectively. We account for foreground reddening using the 
dust maps of lSchlegel et al.l (119981 ) with a value of E(B—V) = 0.04 mag and their prescription 
for converting to HST flight magnitude. Whilst we have encountered much higher reddening 
values in our study of the blue supergiants, the use of only foreground reddening is very likely 
justified in these outer fields. Combining these corrections with the measured magnitudes 
of the tip produces distance moduli of /i = 24.86+o;°J mag for the halo observations and 
H = 24.82+°;^ for with the outer disk field. As an average TRGB distance modulus v 
adopt jj, = 24.84 ± 0.10 mag, which is in basic agreement with that found by IRizzi et al. 
(l2007h using HST WFPC2 data (fi = 24.71 ± 0.04). 



5.3. Discussion 

The two independent methods described in this paper, the FGLR and the TRGB, both 
yield a large distance modulus and agree with each other within statistical uncertainties. The 
FGLR-method stands unaffected by interstellar extinction uncertainties because reddening 
has been determined for each target individually through quantitative spectral analysis, 
which yields accurate stellar parameters and, thus, intrinsic SEDs and colors. In this regard, 
the FGLR work presented here is very similar to the quantita tive spectral analysis of the 



detached eclipsing O-star binary D33 J1013346.2+304439.9 bv iBonanos et all (12009 ). who 
obtained a distance modulus of [i = 24.92 ± 0.12 mag with a reddening of E(B — V) = 0.09 
mag. 

B and A supergiants as well as O stars belong to the young stellar population in the 
dust-obscured disk of M33. An independent and accurate determination of extinction hence 
falls inevitably from their use as distance indicators. The same is also true for other types 
of objects belonging to a population of similar age, such as Long Period Variables and 
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Cepheids. For the former, iPierce et al.1 (120001 ) have carried out a careful photometric study 



of the period-absolute magnitude relationship of LPVs in Per OBI, the LMC, and M33 using 
R and I photometry and a narrow band filter at 8250A (particularly to correct for the varying 
strength of TiO absorption). Assuming a reddening value of E(B — V) = 0.1 mag for M33, 
they obtained a distance modulus of /i=24.85 ± 0.13 mag to the galaxy. While reddening 
was not determined for individual targets, the value assumed is close to the average value of 
0.083 mag found in our study. Thus, it seems consistent that their distance modulus is also 
in agreement with ours. 

The situation with Cepheid distanc e determinations turns out to be more com plex. 



Based on HST WFPC2 photometry from iFreedman et al.1 (ll99ll ) and iLee et al.1 (120021 ). the 



distance modulus was found to be \i = 24.62 ± 0.15 mag and fi = 24.52 ± 0.14 (random) 
± 0.13 (systematic) mag, respectively. Whilst the former study applied a correction for 
metallicity, the latter did not. Both works considered a very high value of reddening, E(B — 
V) = 0.20 mag, which was derived from the difference in the apparent distance moduli 
between V- and /-bands. This reddening value differs substantially from the average value 
of 0.08 mag found for the B and A supergiants in our study. 

Since Cepheids are also young massive stars (though less massive than blue supergiants), 
they are found at similar sites as are blue supergiants, and thus, the reason why their 
reddening is systematically higher is unclear. The direct determination of reddening from 
spectroscopy should be superior to the indirect way of using the difference of apparent 
distance moduli in V- and /-bands. We note that the latter must also depend on metallicity 
(an effect that had not been considered in the pertinent studies) as well as an assumed 
value of reddening for the calibration sample of Cepheids in the LMC. If one assumes that 
the value of E(B — V) = 0.2 mag as adopted for the Cepheids study is too high and uses 
the average value of 0.083 mag from the blue supergiants study instead, then the Cepheid 
distance modulus would increase by 0.37 mag, yielding a distance in excellent agreement with 
our FGLR and TRGB values. This drastic improvement illustrates two things: first, the HST 
observations of Cepheids seem to be compatible with a large distance modulus; second and 
most importantly, the accurate elimination of the reddening uncertainty is absolutely crucial 
for the determination of distances using Cepheids. 



Very recently, IScowcroft et al.1 (120091 ) presented ground-based B- } V-, /-band photome- 
try of Cepheids in the center of M 33 and in the southern spiral arm, some 4 to 5 kpc away 
from the center, with 91 and 28 objects in the two fields, respectively. The use of the We- 
senheit Wyi magnitude, which was assumed to be reddening- free, yielded largely discrepant 
distance moduli in the center and in the spiral arm (24.37 and 24.54 mag, respectively) that 
the authors attributed to a metallicity difference between the regions. They derived a true 
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distance modulus of /i = 24.53 ± 0.11 mag accordingly, very much in agreement with the HST 
photometry Cepheid papers discussed in the previous paragraph. The paper does not pro- 
vide an estimate of E(B — V), but given their distance agreement with the published HST 
Cepheid work, we assume that their reddening value would be close to 0.2 mag. We also 
note that, at least from the basic principle, the use of the Wesenheit magnitude is very sim- 
ilar to the use of difference between apparent distance moduli obtained in V- and /-bands. 
Therefore, we feel that the same discussion on reddening and distance determination in the 
previous paragraphs applies here as well. (This distance modulus is based on a distance 
modulus of 18.40 mag for the LMC, while our FGLR calibration adopted 18.50 mag. Thus, 
it would increase by 0.1 mag if our distance were assumed.) 

It is also impor tant to take int o acco unt that the metallicity difference between the 
two fields studied by lScowcroft et all (120091 ) seems to be significantly smaller than the value 
of 0.566 dex they adopted. This value wa s obtained ba s ed on the assumption of a two- 
component slope as originally proposed by IVilchez et all (119881 ) and then also adopted as 
one possibility by iMagrini et all (j2007bl ). However, as discussed in §4.31 and shown in our 
Fig. [TU both the metallicities from the stars and from H II regions do not support a very 
steep gradient towards the center of M33. Using our result from §4.31 for the abundance 
gradient, the metallicity difference between the two Scowcroft et al. fields is only 0.31 dex. 
With this number, the dependence of the Cepheid distance moduli on logarithmic metallicity 
ch anges would becom e 7 = 0.55 mag dex" 1 , almost twice as high as the value of 0.29 found 
by iMacri et all (120061 ) and that derived by IScowcroft et all (120091 ) . 



There were also a number of other distance determinations methods in the literature 
that involved stars of significantly older populations. The TRGB method, as app lied in this 



work, is a typical example. Our distance modulus agrees very well with that of iKim et al. 
J2002h . who u sed HST WFPC p hotometry in 10 fields and obtained li = 24.81 ±0.13 mag. 



Note also that IKim et all (120021 ) were the first to draw attention to the fact that the shorter 
distance modulus associated with Cepheid distances might be a result of overestimating of 
interstellar reddening. 

In addition, three studie s involving ground-based p hotometry were publish ed later yield- 
i ng qu ite a range in distance: iMcConnachie et all (120041) . \i = 24.50 ± 0.06 mag; iGalleti et al. 
koo4 ). \i = 24.64 ± 0.15 mag; and iTiede etaD J2004T ). /i = 24.69 ± 0.07 mag. It seems 



that systematic effects such as J-band calibration and the algorithm used for TRGB edge 
detection are of importance when comparing these results. We note, however, that two of 
these studies agree with the results acquired with HST photometry within the corresponding 
error margins. 



Other independent investigations using older stellar populations were carried out by lSarajedini et al 
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( 120001 ) and lCiardullo et al.l (120041 ) and confirmed a large distance modulus. The former stud- 
ied Horizontal Branch stars and found /i = 24.84 ± 0.16 mag, whereas the latter used the 
Planetary Nebulae Luminosty Function to obtain fj, = 24.86^oii mag. 



Sarajedini et al.l (120061 ) used HST ACS photometry of RR Lyrae field stars because at 
minimum light, the intrinsic V — I color of RR Lyrae stars is well defined and independent of 
metallicity and period. This allows for an estimate of reddening from the construction of V 
and / light curves. Unfortunately, the reddening determination was somewhat uncertain and 
18 stars in their sample turned out to have negative reddening. The authors found // = 24.67 
± 0.07 mag while excluding those objects, and fi = 24.76 ± 0.08 mag with them included. 
We believe that these 18 stars are simply cases with very small reddening and represent the 
uncertainty of the procedure; thus, their inclusion is justified. This, then, exhibits agreement 
with the HST TRGB results as well as with that derived from the HB and PN stars. 

In summary, we conclude that a large distance modulus as obtained with the FGLR 
and TRGB work presented here is well supported by the value of reddening we found and 
agrees well with other independent work based on low-mass stars of older populations. 



6. Conclusions 

Motivated by the need to achieve precision in extragalactic distance determination 
and to resolve the large discrepancy in the distance modulus to the nearby Triangulum 
Galaxy M33, we undertake an independent investigation applying the flux- weighted gravity- 
luminosity relationship to a sample of 22 blue supergiants at various galactocentric distances 
in M33. With medium resolution spectra from Keck, we carry out a quantitative spec- 
tral analysis and measured stellar parameters such as rotational velocities, stellar gravities, 
effective temperatures, and metallicities. Together with state-of-the-art non-LTE model 
atmospheres spanning an expansive grid in (T e ff, log g) parameter space and photometry 
available in literature, we further derive interstellar reddening, stellar luminosities, radii, and 
masses. The combinations of these parameters allow us to determine a distance modulus of 
24.93 ± 0.11 mag to M33. 

We discuss in detail the stellar properties of our supergiants in relation to their pre- 
dicted evolutionary paths. On the high-luminosity end, our spectroscopic results are nicely 
confirmed by the evolutionary tracks on the H-R diagram. On the low-luminosity end, how- 
ever, we observe one object (similar to the one found by K08) that appears to follow a 
different path and offer close binary evolution as an explanation, which will require further 
confirmation from detailed evolutionary calculations in the future. Furthermore, we note a 
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small systematic discrepancy between spectroscopic and evolutionary masses, however it is 
unlikely to be caused by uncertainties in the distance modulus. Our discussion on metallici- 
ties specifies the significant metallicity gradient as a function of galactocentric distance that 
we find in M33, with large scatter about the regression curve. This is similar to the results 
of several other metal abundance studies of H II regions in literature. 

Our FGLR for M33 follows that for other galaxies very nicely, giving confidence to 
both the reliability of its calibration and the distance modulus thus determined. Using 
HST photometry available in the archive, we also present a TRGB distance modulus that 
nicely converges with our FGLR result within the uncertainties. These distances agree well 
with other published large distance moduli and low reddening values found in literature, 
based on work on detached eclipsing binary, long period variables, horizontal branch stars, 
and planetary nebulae luminosity function. The discrepancy between our results and that 
from Cepheid studies can be explained by differences in interstellar reddening assumed in 
the different work, and we point out that the Cepheid distance moduli will increase if the 
reddening value and metallicity were properly determined and applied. However, we note 
that a final conclusion will only be possible after a more rigorous calibration of the FGLR 
using a large sample of LMC blue supergiants. 

The data presented herein were obtained at the W.M. Keck Observatory, which is oper- 
ated as a scientific partnership among the California Institute of Technology, the University 
of California and the National Aeronautics and Space Administration. The Observatory was 
made possible by the generous financial support of the W.M. Keck Foundation. The authors 
wish to recognize and acknowledge the very significant cultural role and reverence that the 
summit of Mauna Kea has always had within the indigenous Hawaiian community. We are 
most fortunate to have the opportunity to conduct observations from this mountain. 
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Fig. 1. — Comparison of pre- normalized flux spectra for two selected supergiants of late 
(top) and early (bottom) B types reduced by DEEP2 IDL pipeline and by IRAF reduction 
packages, respectively. Notice that the SNR in this blue wavelength range is much better 
for the IRAF-reduced spectra than for the pipeline results. 
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Fig. 2. — The effects of different rotational velocities on the shape of the Mg II (4481) line. 
As the velocity increases, the line grows flatter and flatter, with more pronounced effects in 
the core than in the wings. 
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Fig. 3. — These panels show the Balmer lines fit (from left to right: H 7 , H^, H 3 _ 14 ) for the 
object No. 11. The observed spectrum (thin) is overplotted with the model (thick) with T e ff 
of 9600 K and log g at 1.40ig;$ dex where the errors are represented by dashed lines. 
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Fig. 4. — Balmer line fit curve for object No. 11 in the log g-T c g plane (top) and the log gp- 
T e g plane (bottom). The dashed curves correspond to the maximum fitting errors. 
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Fig. 5. — Fit of the Mg I /II ionization equilibrium along the Balmer line fit curve of object 
No. 5. Three models are overplotted. T cS = 8750K and log g = 1.30 (solid curve), T cS 
= 8300K and log g = 1.05 and T eff = 9250K and log g = 1.50 (both dashed). While the 
Mg II line at 4481A is insensitive to model parameter changes along the fit curve in this 
temperature range and depends only on metallicity, the Mgl lines (vertical bars at the right 
hand figure) show a very strong temperature dependence. A temperature of T e g = 8750±250 
is obtained from this plot. 
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Fig. 6.— Line fit of object No. 1. Top: the Hel lines at 4026A, 447lA, 4713A, and 
5876A along the Balmer line fit curve of object No. 1. Three models are overplotted. T e fj = 
11000K and log g = 1.55 (solid curve), T eff = 10000K and log g = 1.35 and T cS = 12000K 
and log g = 1.70 (both dashed). Bottom: Balmer lines (from left to right: Hp, H 7 , FTj) 
at T c fj = 11000K. The gravities are log g = 1.55 (solid) and 1.45 and 1.65 (both dashed), 
respectively. Note that Hg is affected by stellar wind emission in the line core. 
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Fig. 7. — These panels show the metal line fit for the object No. 5. The observed spectrum 
(thin) is overplotted with the model (thick) with T eff = 8750 K and log g = 1.30 dex at 
increasing metallicities. The best-fit metallicity is log (Z/Z & ) = 0.15 (second row). 
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Fig. 8. — The distribution of reddening E(B — V) for our B and A supergiants; the average 
< E(B — V) > for our sample is 0.083 mag. 
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Fig. 9. — Reddening E(B — V) of M33 H II region s as a function of angular galac- 
tic distance (R^ = 35.40' ; Ide Vaucouleurs et al.lll995l ). The data have been taken from 



Rosolowsky fc Simonl (120081 ). see text. The average < E(B — V) > for this sample is 0.11 



mag (not including the three extreme objects with reddening larger than 0.6 mag). 
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Fig. 10. — The M 33 supergiants of this study in the log g-log T eff plane (top) and in the 
H-R diagram (bottom). Early B types are shown as solid squares, late B and A types as 
solid circles. Evolutionary tracks (iMeynet &: Maederl 120031 1 for different ZAMS masses are 
overplotted. 
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Fig. 11. — The observed mass- luminosity relationship compared with the relationships 
(dashed: late B and A supergiants, solid: early B supergiants) obtained from the evolu- 
tionary tracks of Fig. [TTJ Spectroscopic masses derived from stellar gravity and radius are 
used for the M 33 targets of this study. 
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Fig. 12. — Logarithmic ratio of spectroscopic to evolutionary masses as a function of lumi- 
nosity. The upper diagram contains onl y the M33 objec t s of t his study, whereas the lower 
diagram includes the results by K08 and lUrbaneja et al.l (120081 ) . 
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Fig. 13. — Blue supergiant metallicity [Z] as a function of dimensionless angular galactocen- 
tric distance. Circles: late B and A supergiants; squares: early B supergiants. 
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Fig. 14. — Metallicity of blue supergiants, H II regions and Cepheids as a function of di- 
mensionless angular galactocentric distance. A and B supergiants have the same symbols 
as Fig. [13j Logarithmic ox y gen ab undances of H II regions in units of the solar value as 
published by iMagrini et al.l (]2007al ) are plotted as open squares. Logarithmic neon abun- 
dances of H II re gions norma li zed t o the value for B stars in the solar neighbourhood and 
as obtained from lRubin et al.l (120081 ) are shown as la r ge op en diamonds. The metallicity [Z] 



for beat Cepheids as determined by iBeaulieu et al.l (120061 ) are given as crosses. The solid 
line is the regression for the supergiants only, whereas the dashed lines is the regression for 
all objects. 
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Fig. 15. — Blue supergiant metallicity [Z] as a function of dimensionless angular galactocen- 
tric distance compared to the lo garithmic oxygen abundance s (in units of the solar value) of 
the H II regions investigated by iRosolowsky fc Simonl (120081 ). open circles. 
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Fig. 16. — FGLR fits of the blue supergiants in M33. Top: Apparent de-reddened bolometric 
magnitude vs. flux-weighted gravity. Solid circles are late B and A supergiants and solid 
squares are early B supergiants in M33. The solid line is a linear fit as described in the 
text. Bottom: Absolute bolometric magnitude vs. flux- weighted gravity. In addition to 
the M33 targets, ob 



ects from nine other galaxies investigated in the studies by K08 and 
Urbaneja et al.l (120081 ) are also shown. The solid line is the regression FGLR from K08. The 



dashed curve is the stellar evolution FGLR for Milky Way metallicity. 
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Fig. 17.— Observation of the TRGB in M33. The location of the HST ACS halo fields and 
the disk field are shown in the top figure. The bottom figure shows the two color-magnitude 
diagrams obtained from the halo fields and the disk field. The I magnitude of tip of the red 
giant branch is indicated. 
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Table 1. Target Identification 



No. 


Massey et al. (2006) ID 


Alt. ID 


Reference 


1 


J013351. 20+303224.5 






2 


J013337.09+303521.6 






3 


J013340.47+303503.3 






4 


J013340.84+303822.5 






5 


J013341. 36+303629.6 






6 


J013344.27+304247.2 






7 

8 


J013344.43+303843.9 
J013344.81+303217.8 






9 


J013340.55+303158.7 






10 


J013343.26+303153.8 






1 1 
1 1 


j uioo^tu.ou-t-ouz LVi. 1 


n7^ 
U ( OO 


i 


12 


J013344.66+303631.6 


B215a 


2 


13 


J013229.61+303513.3 






14 


J013315.62+302949.3 






15 


J013351. 56+304005.2 






16 


J013300.23+302323.7 


117 A 


2 


17 


J013339.08+302010.7 






18 


J013300.86+303504.9 


B38, OB 21-108, UIT030 


2, 3,4 


19 


J013342.06+302142.3 


B287, OB 112-41 


2, 3 


20 


J013335.76+310046.9 


B157, UIT136 


2,4 


21 


J013333.72+304719.9 


UIT 122 


4 


22 


J013327.35+310056.4 


UIT 103 


4 



Note. Refe rences : (1) llvanov et al l fll993T) . (2) 

Humphreys fc Sandagel Jl98ph . (3) IMassev et alT Jl995h . (4) 
Massev et all Jl996h 
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Table 2. Spectra, Spectral Types, and Spectroscopic Parameters 



No. 


Sd 


S 1 


v sin i 


C 2 




3 rp4 


- 1 eft 




log g 


\z\ 




log Qf 




Tvnp 




(km s~M 


(km s~ 1 






(K) 




(dex) 




(dex) 




(dex) 




1 


B9 


D 


40 




0.23 


b 


11000±500 


1 


55±0 


14 


-0.07±0.15 


1 


38±0 


06 


2 


A2 


D 


27 




0.15 


a. 


8500±225 


1 


15±0 


15 


0.01±0.12 


1 


43±0 


10 


3 


B9 


D 


53 




0.14 


b 


11500±500 


1 


70±0 


13 


0.08±0.20 


1 


46±0 


06 


4 


A2 


D 


25 




0.09 


a. 


8500±225 


1 


50±0 


20 


0.23±0.15 


1 


78±0 


06 


5 


A2 


D 


22 




0.11 


a 


8750±250 


1 


30±0 


18 


0.10±0.15 


1 


53±0 


13 


6 


B3 


D 


45 




0.14 


c 


16000±1000 


2 


15±0 


20 


0.00±0.10 


1 


33±0 


09 


7 


B9 


D 


30 




0.06 


b 


11000±500 


1 


65±0 


15 


0.12±0.23 


1 


48±0 


08 


8 


B8 


D 


56 




0.22 


b 


12500±500 


1 


90±0 


12 


-0.01±0.12 


1 


51±0 


05 


9 


B3 


D 


58 




0.23 


c 


16000±1000 


2 


15±0 


10 


0.00±0.10 


1 


33±0 


09 


10 


AO 


D 


42 




0.23 


a,b 


10000±500 


1 


35±0 


15 


-0.02±0.20 


1 


35±0 


06 


11 


AO 


E 


35 


30 


0.54 


a,b 


9600±340 


1 


38±0 


11 


-0.50±0.18 


1 


45±0 


05 


12 


AO 


E 


50 


30 


0.09 


a,b 


9700±270 


1 


ui/XU 


J. _ 


-0.05±0.21 


1 


O I _ 1_U 


D7 


13 


AO 


E 


42 


20 


0.76 


a,b 


9550±320 


1 


41±0 


12 


-0.48±0.19 


1 


49±0 


06 


14 


B9 


E 


35 


25 


0.37 


a,b 


10000±410 


1 


32±0 


16 


-0.37±0.21 


1 


32±0 


09 


15 


B9 


E 


45 


30 


0.01 


a,b 


11500±500 


1 


63±0 


16 


0.03±0.22 


1 


39±0 


08 


16 


AO 


E 


45 


30 


0.57 


a,b 


9600±328 


1 


10±0 


12 


-0.52±0.19 


1 


17±0 


06 


17 


B3 


E 


50 


40 


0.59 


c 


17000±1000 


2 


00±0 


15 


-0.17±0.20 


1 


08±0 


05 


18 5 


Bl 


I 


60 




0.46 


c 


22000±500 


2 


60±0 


10 


-0.22±0.15 


1 


23±0 


14 


19 5 


BO 


E 


65 




0.54 


c 


26000±500 


2 


79±0 


10 


0.00±0.15 


1 


13±0 


13 


20 5 


BO. 5 


E 


70 




0.74 


c 


24000±500 


2 


60±0 


10 


-0.30±0.15 


1 


08±0 


13 


21 5 


B0.5 


E 


70 




0.35 


c 


24000±500 


2 


77±0 


10 


-0.22±0.15 


1 


25±0 


13 


22 5 


B0.7 


E 


50 




0.79 


c 


23700±500 


2 


65±0 


10 


-0.52±0.15 


1 


15±0 


13 



Spectrograph: D - DEIMOS, E - ESI, I - ISIS 
2 Macroturbulence velocity 



3 R25 = 35.40', the radius of M33 in the plane of the galaxy; this has been corrected for a 
position angle of 22° and an inclination angle of 54° 

4 T efT -method. a: Mg i/ll, b: He I, c: Si Il/lIl/lV 

5 Parameters adopted from lUrbaneja et al.l (l2005bl ). with the galactocentric distances cor- 
rected for the normalization radius and the distance to the galaxy used in this paper 
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Table 3. Photometric Properties and Bolometric Corrections 



No. 




my 1 




B - V 1 


(B - V) 


E(B - 


V) 




My 2 




BC 






(mag) 




(mag) 


(mag) 






(mag 






(mag" 




(mag) 


1 


17 


194±0 


004 


0.034±0.004 


-0.046±0.021 





080±0 


021 


-7 


984±0 


065 


-0.481±0.082 


2 


17 


937±0 


004 


0.062±0.004 


0.042±0.032 





020±0 


032 


-7 


055±0 


099 


0.005±0.030 


3 


17 


997±0 


004 


0.025±0.004 


-0.067±0.020 





092±0 


020 


-7 


218±0 


062 


-0.555±0.085 


4 


18 


235±0 


003 


0.172±0.004 


0.012±0.025 





160±0 


025 


-7 


191±0 


078 


0.025±0.045 


5 


17 


913±0 


004 


0.106±0.004 


0.019±0.028 





087±0 


028 


-7 


284±0 


087 


-0.032±0.032 


6 


18 


209±0 


004 


0.011±0.004 


-0.131±0 


022 





142±0 


022 


-7 


161±0 


068 


-1.343±0.139 


7 


17 


941±0 


003 


0.054±0.003 


-0.055±0 


020 





109±0 


020 


-7 


327±0 


062 


-0.458±0.085 


8 


18 


067±0 


004 


0.025±0.004 


-0.093±0 


013 





118±0 


014 


-7 


229±0 


044 


-0.740±0.087 


9 


18 


047±0 


004 


-0.039±0.004 


-0.131±0 


010 





090±0 


011 


-7 


162±0 


034 


-1.343±0.070 


10 


17 


221±0 


004 


0.046±0.004 


-0.011±0 


023 





057±0 


023 


-7 


886±0 


071 


-0.301±0.081 


11 


17 


856±0 


004 


0.044±0.004 


-0.016±0 


020 





060±0 


020 


-7 


260±0 


032 


-0.240±0.050 


12 


16 


904±0 


004 


0.031±0.004 


-0.024±0 


020 





055±0 


020 


-8 


197±0 


032 


-0.240±0.050 


13 


17 


623±0 


004 


0.108±0.004 


-0.020±0 


009 





130±0 


010 


-7 


710±0 


031 


-0.230±0.040 


14 


17 


120±0 


004 


0.061±0.004 


-0.019±0 


030 





080±0 


030 


-8 


058±0 


093 


-0.330±0.050 


15 


17 


102±0 


003 


0.039±0.003 


-0.081±0 


030 





120±0 


030 


-8 


200±0 


093 


-0.570±0.060 


16 


16 


440±0 


004 


0.142±0.004 


-0.012±0 


020 





130±0 


030 


-8 


893±0 


093 


-0.300±0.020 


17 


17 


209±0 


005 


-0.075±0.005 


-0.150±0 


019 





075±0 


020 


-7 


954±0 


032 


-1.490±0.140 


18 


17 


322±0 


004 


-0.118±0.004 


-0.170±0 


010 





052±0 


011 


-7 


769±0 


034 


-2.110±0.083 


19 


17 


989±0 


004 


-0.179±0.004 


-0.190±0 


010 





012±0 


011 


-6 


975±0 


034 


-2.470±0.084 


20 


17 


898±0 


005 


-0.144±0.005 


-0.170±0 


010 





026±0 


011 


-7 


113±0 


034 


-2.340±0.084 


21 


17 


754±0 


004 


-0.063±0.004 


-0.180±0 


005 





117±0 


006 


-7 


539±0 


019 


-2.340±0.084 


22 


18 


035±0 


005 


-0.154±0.005 


-0.180±0 


005 





026±0 


007 


-6 


976±0 


022 


-2.320±0.084 



1 Adopted from lMassey et al.l (120061 ) 

2 Computed with distance modulus of 24.93 ± 0.11 mag (this work) 
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Table 4. Bolometric Magnitudes, Luminosities, and Stellar Masses 



1NO. 




M bol l 


10] 


1 (L/L Q ) 


it 


A/f 

-'"spec 




M evol 








(mag) 






(dex) 










(M Q ) 




1 


-8. 


.51±0. 


.11 


5. 


30±0. 


.04 


121±13 


19.0±8. 


,8 


20. 


2±0. 


.7 


2 


-7. 


.05±0. 


10 


4. 


73±0. 


.04 


106±8 


5.8±3. 


.2 


13. 


1±0. 


.4 


3 


-7. 


.77±0. 


.11 


5. 


01±0. 


.04 


81±8 


11.9±5. 


.1 


16. 


2±0. 


.5 


4 


-7. 


17±0. 


09 


4. 


77±0. 


.04 


112±8 


14.3±8. 


.8 


13. 


6±0. 


.3 


5 


-7. 


,32±0. 


09 


4. 


83±0. 


.04 


113±8 


9.3±5, 


.1 


14. 


2±0. 


.4 


6 


-8. 


,50±0. 


16 


5. 


30±0. 


,06 


58±9 


17.5±12. 


.5 


20, 


9±1. 


.0 


7 


-7. 


.79±0, 


.11 


5. 


02±0. 


.04 


89±10 


12.8±6. 


.3 


16. 


3±0. 


.5 


8 


-7. 


,97±0. 


.11 


5. 


09±0. 


.04 


75±7 


16.1±6. 


3 


17. 


2±0. 


.5 


9 


-8. 


k /i _i_n 
,o4±U. 


no 


5. 


on in 

oz±U. 


no 


58±4 


17.5±5. 


,4 


20. 


9±0. 


.5 


10 


-8. 


19±0. 


.11 


5. 


18±0. 


.04 


129±15 


13.6±6. 


.9 


18. 


4±0. 


.6 


11 


-7 


.51±0. 


.08 


4. 


90±0 


.03 


103±9 


9.2±3. 


.2 


15. 


0±0. 


.3 


12 


-8. 


.44±0. 


.08 


5. 


28±0. 


.03 


154±11 


18.0±6. 


.5 


20. 


0±0. 


.5 


13 


-7. 


,94±0. 


.05 


5. 


08±0 


,02 


127±9 


15.0±5. 


.4 


17. 


1±0. 


.3 


14 


-8. 


.39±0. 


.11 


5. 


26±0. 


.04 


142±14 


15.2±7. 


,8 


19. 


7±0. 


.7 


15 


-8 


,77±0. 


.11 


5. 


41±0. 


,04 


128±13 


25.4±13. 


.2 


22. 


4±0. 


.8 


16 


-9 


19±0. 


.07 


5. 


58±0. 


.03 


223±17 


22.8±8. 


,4 


26. 


0±0. 


.6 


17 


-9. 


.44±0. 


.15 


5. 


68±0. 


.06 


80±12 


23.1±12 


.7 


29. 


3±1 


.7 


18 


-9. 


.88±0. 


09 


5. 


85±0. 


.04 


58±4 


48.9±14. 


.6 


35 


2±1. 


.4 


19 


-9 


,45±0. 


09 


5, 


68±0. 


.04 


34±2 


26.1±7. 


.6 


29. 


3±1 


.0 


20 


-9 


.45±0. 


09 


5. 


68±0. 


.04 


40±2 


23.3±6. 


.9 


29. 


,4±1. 


.0 


21 


-9. 


.88±0. 


09 


5. 


85±0. 


.04 


49±3 


51.1±15 


.0 


35. 


2±1. 


.3 


22 


-9 


.30±0. 


09 


5. 


62±0. 


.04 


38±2 


23.8±7. 


.0 


27. 


7±0. 


.9 



Computed with distance modulus of 24.93 ± 0.11 mag (this 
work) 



